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Abstract 

We analyze some solvable models of a quantum mechanical system in interac- 
tion with a reservoir when the initial state is not factorized. We apply Nakajima- 
Zwanzig's projection method by choosing a reference state of the reservoir endowed 
with the mixing property. In van Hove's limit, the dynamics is described in terms of 
a master equation. We observe that Markovianity becomes a valid approximation 
for timescales that depend both on the form factors of the interaction and on the 
observables of the reservoir that can be measured. 
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1 Introduction 



The dissipative dynamics of a small quantum system weakly coupled to a large 
reservoir is described in terms of a master equation [1,2,3,4]. In the standard 
approach to this problem, one usually takes for granted that there are no initial 
correlations between the system and the reservoir. In the preceding article [5], 
hereafter referred to as Article I, we reconsidered this hypothesis in the frame- 
work of Nakajima-Zwanzig's projection method [2,4,6,7] and proved that, in 
order to get a consistent description, the reference state of the reservoir should 
be endowed with the mixing property. In such the initial correlations 

disappear in the Markovian (van Hove) limit and the system behaves as if 
it started from a factorized initial condition. Interestingly, one arrives at the 
same conclusions also for uncorrelated initial conditions. The mixing property 
is therefore crucial, and a "wrong" choice of the reservoir state provokes the 
appearance of secular terms. 

In this article, we shall focus on the hypotheses that are necessary for the 
derivation of the theorem proved in Article I [5] . These will be scrutinized in 
terms of two exactly solvable models, in which an oscillator is coupled to a 
bosonic reservoir. This will enable us to describe the onset to Markovianity 
and the timescales at which Markovianity becomes a valid approximation. 

This article is organized as follows. We introduce notation and summarize 
previous results in Sec. 2. The first exactly solvable model is introduced in 
Sec. 3 and solved in Sees. 4-6. The second model is briefly discussed in Sec. 
7. Section 8 is devoted to a discussion and some concluding remarks. Two 
Appendices contain the details of the derivations. 



2 Summary of Previous Results 



2. 1 Notation 



We start by briefly summarizing the main ideas of Article I [5] and introduce 
notation. Let the total system consist of a "large" reservoir B and a "small" 
(sub)system S, so that the total Hilbert space can be expressed as the tensor 
product of the Hilbert spaces of the reservoir H B and of the system Hs, 

H tot = H s ®H B . (2.1) 

The Hamiltonian and the corresponding Liouvillian of the total system read 

H = H + XH SB = H S + H B + A#sb, (2.2) 
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C = C + A£ SB = C s + C B + A£ S b, (2.3) 
respectively, where A is the coupling constant. Clearly, 

[H s ,H B ] = 0, [£ S ,£ B ]=0. (2.4) 

We assume that the system Hamiltonian H s admits a pure point spectrum, 
and the system Liouvillian £ s is resolved in terms of its eigenprojections Q m , 

C-S i ^ ' ^mQrni ^ ' Qrn lj QmQn ^rnnQrri' (^-5) 



2.2 Nakajima-Zwanzig 's Projection Method 



Let p(t) be the density operator of the total system at time t, which has 
evolved from the initial state p 

p(t) = e a Po (2.6) 
and is the solution of the von Neumann equation 

j t p{t)=Cp{t), p(0)= Po . (2.7) 

We are interested in the reduced dynamics of system S, which is described by 
the density operator of S, 

p s (t) =tr B p(t). (2.8) 

In order to derive a master equation for ps{t), Nakajima-Zwanzig' s procedure 
makes use of the projection operators [2,4,6,7] 

Vp = tr B {p} <g> n B = a® Q B , Q=l-V, (2.9) 

where fl B is a certain reference state of the reservoir. Due to normalization 
tr B I^b — 1, it follows that V 2 = V and Q 2 = Q. In particular, 

v P (t) = Ps (t)®n B , Qp(t) = p(t)-p s (t)®n B , (2.io) 

where we used the definition (2.8). 

In the standard derivation of a master equation, the initial state of the total 
system, p , is taken to be the tensor product of a system initial state ps an d 
a reservoir state pb, 

Po = Ps®Pb- (2.11) 
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This is an uncorrelated initial state. The reservoir is assumed to be at equi- 
librium (with respect to the reservoir free evolution £ B ) 

CbPb = 0, (2.12) 

and in most applications p B — ZZ x e~® H,R is a thermal state at the inverse 
temperature f3 = (A; B T) _1 with the normalization constant Zp. Then, the 
reservoir state p B in the uncorrelated initial state (2.11) is usually taken as 
the reference state fi B . 

When the assumption of a factorized initial state is not justified, however, an 
ambiguity arises regarding the choice of the reference state f2 B . Indeed, if 

Po = Ps ® Pb + 5p , (2.13) 

where 

Ps = tr B po, p B = tr s Po, (2.14) 

and the term Spo represents the correlation between system S and reservoir 
B, the relation between p B and fi B is by no means obvious. We discussed this 
point in Article I [5] and proved the following theorem. 



2.3 Theorem 



Given a correlated initial state po, if 

(i) is the unique simple eigenvalue of the reservoir Liouvillian £ B corre- 
sponding to the eigenvector f2 B and the remaining part of the spectrum of 
Cb is absolutely continuous (strictly speaking, the spectrum of £ B can be 
defined only once the sector has been specified: in our case, the relevant 
sector is that containing the state f2 B ); 

(ii) the initial (correlated) state of the total system is given in the form 

p = A(l s <g) n B ) = ]T L,(l s <g> n B )L\, (2.15) 

i 

where A is a bounded superoperator (i.e., Lj's are bounded operators) 
satisfying the normalization condition trp = 1, namely, the initial state 
po belongs to the sector specified by Is <8> ^ B , 

then van Hove's "A 2 t" limit [1,8,9] of the P-projected density operator in the 
interaction picture, 

Pl (r) = limp| A) (r) = ^e'^V p{r / X% (2.16) 

A — >0 A — >U 
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is the solution of 



Pi(r)=V Po + J dr'JCp^r') 



(2.17) 



o 



with 



Q 



£sBQm'P, 



(2.18) 



rn 



C + iu m - 0+ 



or equivalently, 



d_ 




/Cpi(r), 



pi(0) = Vpo = tr B {p } <S> ^b- 



(2.19) 



That is, even if the initial state p is n °t in a factorized form, all correlations 
disappear in van Hove's limit and system S behaves as if the total system 
started from the factorized initial state in (2.19) with the reservoir state Qb- 

In addition, we showed that 



which makes the dynamics consistent, for no spurious term will develop in the 
master equation and no correlations can appear at later times: not only the 
initial state, but also the state at any moment t is factorized in van Hove's 
limit. This supports the validity of the assumption of the factorized state, that 
is frequently applied in literature in order to derive a master equation [2,3,4]. 
The state of system S evolves according to the master equation (2.19), while 
the reservoir B remains in the state 

It is important to note that, in van Hove's limit, the reservoir state imme- 
diately relaxes into Qb, which is the eigenstate of the reservoir Liouvillian 
£b belonging to its unique simple eigenvalue 0, and the spectral properties 
required in hypothesis (i) imply that it is a mixing state. The right choice for 
the reference state of the projection (2.10) is this mixing state Qb, and such 
a projection is nothing but the eigenprojection of the reservoir Liouvillian jCb 
belonging to the simple eigenvalue 0. This is the criterion for the reference 
state, that covers both equilibrium states and nonequilibrium steady states. 
Furthermore, as clarified in Article I [5], the interaction between system S and 
reservoir B is not essential to the factorization or the mixing; the total system 
is factorized and the reservoir relaxes into the mixing state through its own 
free evolution. 

The purpose of the present article is to scrutinize these issues in some explicit 
examples. In particular, we shall focus on: (a) the disappearance of the initial 
correlation, (b) the factorization of the total system, and (c) the relaxation of 
the reservoir into the mixing state, in van Hove's limit. This will also enable 
us to discuss the relevant timescales for the factorization and the mixing. 



lim Qp(r/A 2 ) = 0, 



(2.20) 
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3 An Exactly Solvable Model 



Let us corroborate the above general arguments by scrutinizing an exactly 
solvable model. We consider an oscillator a coupled to a reservoir b^, whose 
Hamiltonian is given by (2.2) with 

oo oo 

H s = w s a f a, H B = J duublb u , H SB = i j duj {gla ] b u - g u abl), (3.1) 
o o 

where a (a*) and b w (b^) are annihilation (creation) operators satisfying the 
canonical commutation relations 

[a,at] = l, [b u ,bl] = 6(u-u'), (3.2) 

and g w is the form factor of the interaction. Even though system S has an 
infinite number of levels, and does not fulfill the conditions of the main theorem 
proved in Article I [5], the following explicit calculation will show that all the 
conclusions are still valid and therefore the theorem has a wider applicability. 

The above model is exactly solvable [10,11,12]. Indeed, the Heisenberg equa- 
tions of motion for a(t) = e tHt ae~' lHt and b w (t) = e^b^e' 1111 read 

oo 

a(t) = —iu) S a(t) + \ J duj g*b u (t), (3.3a) 
o 

Kit) = -iub u (t) - Xg u a(t), (3.3b) 

and by integrating the second equation and inserting it into the first, one 
obtains an integro-differential equation for a(t), 

t 

a{t) = -iu s a(t) -X 2 J dt' Kit - t')a{t') + XB(t), (3.4) 

o 

with 

oo oo 

Kit) = Jdw \gj 2 e^\ Bit) = J dwg^X, (3.5) 



which is solved via Laplace transform to yield 

t 

ait) = G(t)a + A J ' dt' G(t - t')B(t'), (3.6a) 
o 

t 

b u (t) = e'^X -Xj dt' e-^-^&X*'), (3.6b) 
o 
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where 

G(t) = [ ^ 6 ——, , K{s) = Jdcu-^-, (3.7) 

C B being the Bromwich path on the complex s-plane. Note that G(0 + ) = 1 
and G(0+) = -iu s . 



4 A Correlated Initial State 



Any physical preparation of a quantum state is based on concrete physical 
procedures that cannot be controlled with complete accuracy. The real initial 
state is therefore unknown to some extent and in general has certainly some 
correlations built in. As an example of a correlated initial state, that has the 
advantage of being solvable, we take 

Po = ^^(a s ®p w )e^ a , (4.1) 

with any positive operator a s of system S and a reservoir state 

Pw = ^e-^\ (4.2) 

where the summations over the reservoir modes u are implicit (and so hence- 
forth as long as no confusion can arise): 

oo oo oo 

?b = JdwCbu, tfWb = Jduj J dw'blW^. (4.3) 



W^c/ is Hermitian (W<^y = and consists of Wf®,, that is proportional 

to 5{uj — u)'), and the remaining square integrable part VV wt y, 



W<*s = n£, + W^, = W{u)5{u - J). (4.4) 

The states po and pyy are normalized with the normalization constants Zq and 
Zyy, and ^ is the relevant parameter to the initial correlation between system 
S and reservoir B. 

F° r £,u — 0, the state (4.1) is obviously factorized, while it becomes a tightly 
correlated state for any £ w ^ 0, with correlations proportional to as will 
be shown later in (5.7). Actually, the operator e at ^ tfe appearing in the ini- 
tial state (4.1) generates a correlation between S and B: it changes the n- 
particle states Wrji ■ ■ ■ 6^ n |vac) of the reservoir into e^^^Wrji ■ ■ • 6^ n |vac) = 
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(tfrji + a^^rji) ■ ■ ■ (b^rjn + <rf^r) n ) |vac), so that S and B are entangled for any 
nonvanishing value of £ w . It is also possible to explicitly compute the correla- 
tion functions in the initial state (4.1): see the generating functional (5.4) and 
the correlation function (5.7) below. The choice of this particular form for the 
initial state p is mainly due to the fact that it allows us to solve the dynamics 
of the total system exactly and to discuss the correlation between system S 
and reservoir B. One can think of the correlations in (4.1) as engendered by 
a linear interaction of the form H prep oc aJ^b + h.c. in a rotating- wave-like 
approximation. 

As shown in Appendix A of Article I [5], the (normalized) reservoir state 

n B = -L e -^ 0) » (4.5) 

is mixing with respect to the reservoir dynamics driven by the Hamiltonian 
H-q in (3.1), and the initial state p m (4-1) belongs to the sector specified by 
l s <E> in the sense of (2.15). Indeed, p is the state perturbed from l s ® H B 
by a local operator L, 



where 




is a local perturbation such that 

p w = Lb^b^b, (4.8) 



T denoting the anti-chronologically ordered product and 

oo oo 

6t e -^ (0) We^ (0) 6 = J duJ J duJ ' ble'^W^e^'K, (4.9) 

o o 

Even though the initial state po does not satisfy the hypotheses of the theorem 
proved in Article I [5], the following analysis extends the general results valid 
for a bounded perturbation. 

Note that the reservoir Gaussian state pw in (4-2) is fully characterized by 
the two-point function 

ACu/ = (&1A)pw = tr B {&lAPw} 
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and, as shown in Appendix A in Article I [5], it is also composed of two parts 
like in (4.4), 

AL- =M ( 2,+N^,. (4.11) 
The first term is the two-point function in the mixing state 

^ = <blbJ) nB =N(u)6(u-cS), N{u) = eW * ) _ 1 , (4.12) 

which is the Bose distribution function when W(u>) = f3u, while the second 
one is a local function representing the effect of the local perturbation Lb in 
(4.8). 



5 Dynamics of the Total System 

Since we are interested in the correlation between system S and reservoir B, 
we need to look at the state of the total system, p(t). In order to treat the 
reservoir degrees of freedom rigorously, we should restrict ourselves to reservoir 
observables whose expectation values are finite and discuss the state of the 
total system through a set of such expectation values. The relevant quantity 
for our discussion is therefore a characteristic functional of the state p(t), e.g. 

G[J a , J* a , J b , J\- 1] = tr{e J >e J ^e- 6t ^e- atJ v(0}, (5.1) 

where j\b = / °° duo J£ u b u , which is the generating functional of the expec- 
tation values of any anti- normally ordered products of a, a\ b u , and &t and 
characterizes the state of the total system, p(t). It is important to note that 
we are not interested in infinitely extended objects, such as the Hamiltonian 
of the reservoir Hb, since their expectation values are infinite: our targets are 
locally distributed objects. Such a formalization of the problem is reasonable, 
since we cannot observe infinitely extended objects in practice, and this is 
nothing but the starting point of the C*-algebraic approach to quantum sta- 
tistical mechanics [13]. In the characteristic functional (5.1), the bandwidth 
of J b:U! represents the locality of the observables. 

Let us begin with the characteristic functional of the initial state po in (4.1), 
G [J a , J* a , J b , f b ] = g[J a , J* a , J b , f b - 0] = tr{e J >e J t b e- bU »e- aU y }, (5.2) 

which, in the coherent-state representation (Q-representation [4]) 
a\a) = a\a), (a\a') = e -M 2 /2-M 2 /24W f fa = ^ 

J 7T 

is evaluated as 
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Qo[Ja, Ja, Jb, Jb) 

= L e -JU b /^(a| < Ts|a)e J i a - aV -<e 6t « a -^e^ t+J ^) 

Zq J 71 



W 



= J_ e -Jl(l+^)J b f ^ ^ as ^ e a*^a e (J:+jlAfOa e -a*(J a +^MJ b ) 
Zq J TT 

= e- J ^ 1+M ^g s (J a + £WJ b , J* a + 4m, (5-4) 

where 

Gs(J a , J* a ) = Qo[Ja, J* a , 0, 0] = tr s {e J >e- at J Vs} (5.5) 
is the characteristic function of the initial state of system S and 

Ps = tr B po- (5.6) 

One can see from this characteristic functional how the parameter C, w embodies 
the initial correlation. For example, 



d 5 



j j* j, rt_n 

, OO 

/ dJ GAU ^j Ja d jP = J do/ &AU (a« f )ps • 

n a a J„,J*=0 n 

(5-7) 



Ja,J;=o o 



Let us discuss the evolution of the state of the total system 

pit) = e- iHt p e iHt . (5.8) 

The characteristic functional (5.1) of the state p(t) is easily computed in the 
Heisenberg picture 

g\j a , r a , J b , 4; t] = tr{ e ^w e J 2^ e - 6t w j » e - ot w j - A>0 } 

= tT{e Ji(t)a e Jt(t)b e -btj b (t) e -aU a (t) po} 

=g [j a (t),r a (t),j b (t),4(t)\, (5.9) 

where J a (t) and J b (t) are functionals of J a and Jb, defined via a)(t)J a + 
b\t) J b = a)J a (t) + bU b (t). Note that the solutions (3.6) for a(t) and b u (t) 
are linear in a and b w , but do not contain a* or b^. The characteristic func- 
tional of the initial state p is given in (5.4) and Eq. (5.9) is further reduced 
to 

g[j a , j* a , j b , Ji t\ = e - j i^ i+ ^ j ^g s (j a (t) + emit), jm + Jimo- 

(5.10) 
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We thus obtain the exact characteristic functional of the state of the total 
system, p(t), 

g[j a , r a , j b , 4; t] = e- JtA u j g s (h\t)j, j^h(t)), (5.11) 

where 



jU(t)j = (j; 4) [ A ^\ ( J ; ] , u m = (* :(t) *t (t) 

V J Aba(t) A bb {t) I \ T, I 



V 



jj' h ^ t)J= ( h ^^)[ J J b 



(5.12) 
with 



* t 



A aa {t) = X 2 J dt' J dt"G(t - t')$ gg (t',t")G*(t - t"), (5.13a) 







J j b A bb (t)J b = <S> JbJb (t, t) - 2X 2 Re J dt' (K Jb9 *G)(t- t')$ gJb {t' , t) 



t t 

+ A 4 | dt' J dt"(K Jbg *G)(t-t')$ gg (t'X)(G**K* Jbg )(t-t"), (5.13b) 







f b A ba (t) = X J dt'<S> Jbg (t,t')G*(t-t') 


t t 

-X s J dt' J dt" (K Jb9 * G)(t - t')$ gg (t',t")G*(t - t") (5.13c) 


and 

h a (t) = G(t) + X(G * K g(m )(t), (5.13d) 

Jlhit) = K Mm (t) - X(K Jb9 * G)(t) - X\K Jbg * G * K g(m )(t). (5.13e) 

We have introduced 

00 00 00 

K fg (t) = J dufte-^gv, $ fg (t,t') = JdooJ do/ 1 + UU^'g^, 
00 

(5.14) 

where l UU) r = 6(00 — 00'), and the convolution 

t 

(F*G){t) = j dt' F(t-t')G{t'). (5.15) 
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The characteristic functional of the total system (5.11) is exact and valid for 
any time t. 

The functions X(G * Kg(Nt))(t) in h a (t) and Kj b ^^(t) in hb(t) describe how 
the initial correlation propagates, while Ab a (t) and X(K Jbg * G)(t) in hb{t) 
describe the correlation established through the interaction between system S 
and reservoir B. System S forgets its initial state through the decay of G(t) 
and approaches an equilibrium state via the action of A aa (t), while Abb(t) 
governs the relaxation of reservoir B into its equilibrium, i.e. the mixing state 
f^B, as explained in the following. 



6 The van Hove Limit of the Characteristic Functional and Dis- 



We are now in a position to discuss the van Hove limit of the evolution of the 
total system and demonstrate the validity of the general theorem proved in 
Article I [5]: (a) the disappearance of the initial correlation, (b) the factoriza- 
tion of the total system, and (c) the relaxation into the mixing state, in van 
Hove's limit. 

In order to discuss van Hove's limit, let us remove the (rapid) oscillation of 
system S. That is, let us look at the characteristic functional of the density 
operator e lHst p(t)e~ lHst in the scaled time r = A 2 t, 



Gi x) [J a , J* a) J b , 4; r] = Q[J a e-^/ x \ jy—' x \ J b , 4; r/A 2 ]. (6.1) 



Then, the van Hove limits of the constituent functions (5.13) (Appendix A), 



cussion 



hm Aa(r/A 2 ) = [1 + N(u; s )](l - e~ r ^), 

A — >\J 



(6.2a) 



oo 




(6.2b) 



o 



\ime luJsr/X2 A ab (T/X 2 )Jb = lim J^A^r/A 2 )^ 3 ^ 2 = 0, 

A — ^0 A — ^0 

lim h a {T/\ 2 )e lu)ST ' X2 = e -r(^W2 e -iA(o;s)r ) / \ 2 ) 

A — >0 X — >0 



(6.2c) 
0, (6.2d) 



lead us to the van Hove limit of the characteristic functional (5.11), 
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&[J a , J a *, J 6 , 4; r] = lim ^ (A) [Ja, J a *, J 6 , 4; r] 

A — 

= e -^V a [l+W(^)](l-e- r ("s)-) e -/ o ° c ^J^[l+7VH]J^ 

x ^ s (^e" r( " s)r/2 e tA(ws)r , j'e-r^W^-^^)'"), 

(6.3) 

where 

r( W ) = 2vr|,J 2 , A( W ) = P J^iVI. (6.4) 

o 

It is clear from (6.3) that (a) the initial correlation (or, equivalently, £ w ) dis- 
appears and (b) the state of the total system is factorized at all times in van 
Hove's limit. Furthermore, (c) the local perturbation in the initial state p , 
i.e. L in (4.6) (especially, the contribution of VV^y, which appears in the char- 
acteristic functional through M UKO i\ decays out and the reservoir relaxes into 
the mixing state Qb given in (4.5). The dynamics of the system in van Hove's 
limit is exactly the same as that derived from the uncorrelated initial state 
tr{p } <8> with the mixing state Qb and it is actually possible to show that 
the density operator pi(r) characterized by the characteristic functional (6.3) 
obeys the master equation 



^:Pi( T ) = -i[A(^ s )ay Pi(r)} 

- i[l + iV(^s)]r(cJs)[a t api(r) + Pl (r)a^a - 2ap I (r)a t ] 

- ^N(u s )T(uj s )[aa^pi(T) + pi(r)aa f - 2a f pi(r) a ]- (6.5) 

This is nothing but the familiar master equation derived from the factorized 
initial condition with the reservoir in the thermal equilibrium state at a finite 
temperature, p ~ ps ® e _/3i?B [3,4]. 

These points corroborate the theorem in Article I, suggesting that the mixing 
state fie, which is contained in the initial state p , should be selected as the 
reference state of Nakajima-Zwanzig's projection V. Note that the character- 
istic functional in van Hove's limit, Eq. (6.3), approaches 

ft[j a , r a , J b , 4; r] ^ e -w+m<*)] e -sr<toj;,jwMv>, u , (6 . 6 ) 

which means that the equilibrium state (in van Hove's limit) is 

p eq = -Le-^ ata ® n B , Z^ = 1 - e~ w ^\ (6.7) 

^cq 

i.e., system S relaxes into the equilibrium state with the same structure as 
that of the mixing state f^- 
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As discussed in Article I [5], the state of the total system is factorized through 
its free evolution and the interaction between system S and reservoir B is not 
essential, which is also confirmed by the present exact solution. In the absence 
of the interaction, the exact characteristic functional (5.11) reads 

g[j a , r a , j b , 4; t] = e -»w*.«>0 s (,7 ae *-B' + K* Mm (t), r a e-^ + K Mm (t)) 

(6.8) 

and approaches 

g[J a ,.r a ,J b ,4;t] ^ e-r^V+^^g^jy^^J^e-^) (6.9) 

by Riemann-Lebesgue's lemma [see the discussion below and Eq. (A. lb)]. The 
state is thus factorized into 

p(t) ^ Ps (t)®n B , p s (t) = e Cst tT BPo (6.10) 

through the free evolution, which confirms the second part of the theorem in 
Article I [5]. 

The timescales of the factorization and the relaxation into the mixing state 
are clear from (6.8): the former is governed by the function 

OO , OO v 

Kj bim (t) = jdujljduj' J 6 yV^£y 1 e-^ (6.11) 

0^0 ' 

contained in jlhb(t) in (5.13e), and the latter by the leading term of Abb(t) in 
(5.13b), 



OO / OO 



*WM) = *S?j 6 (0) 

+ 2 Re J dw ( J dw J£ j u +u /2tf(w+w/2)(w-u//2)Jb,u>-u>/2 | e' luJt , 

(6.12) 

where 

OO 

*g>(f) = J du>f:[l + N(u>)]g u e-**. (6.13) 
o 

The timescales of the decay of these functions are determined by the band- 
widths of their Fourier transforms, 



2tt J dw' Jtjf^& (6.14) 
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for the former, and 

oo 

®J b J b (u) = 2vr J dw JZ fi ^ /2 tf(Q^/2)(B- u /2)Jbfi>-u/2 (6-15) 

w/2 

for the latter. Therefore, besides the spread of the initial correlation and 
of the perturbation from the mixing state (Af wtJ >), the size of the relevant 
reservoir observables (Jb, u ) influences the timescales of the factorization and 
the mixing. In the weak-coupling regime they are very short compared with 
the timescale of the dissipative dynamics of system S, which is of order 1/A 2 
in the original time t. 

It is interesting to discuss what happens from a physical point of view. The 
initial correlations and the local perturbations propagate outwards from the 
region of interest (defined by the "size" of the relevant local observables) and 
never come back. What remains is the "unperturbed" state, that is the mixing 
state f2e and is the stable "ground state" within the sector it specifies. The 
relaxation time of such a process is the time necessary for the disturbance 
to pass through the range of the interaction, that of the initial correlation, 
and the extension of the observable. It should be noted that, when we work 
in the interaction picture e lHot p(t)e~ tHot , instead of e lHst p(t)e~ iHst considered 
above in (6.1), we should duly take into account the time dependence of the 
observables in such a picture, X(t) = e tHot Xe~ lHot , otherwise mixing is not 
observed. 



7 A Solvable Model with Counter-Rotating Interaction 

Let us look at another solvable example: the same model as the previous one 
(3.1) but with a different interaction Hamiltonian 

oo 

H S b = i(a + a f ) J du (g*b u - (7.1) 
o 

containing counter-rotating terms. This model is also exactly solvable [10,12,14,15,16]. 
Let us only briefly sketch the main results. More details are given in Appendix 
B. The exact characteristic functional of the state of the total system, p(t), 
reads 

Q[J a , J* a , J b , 4; t] = e -J^(t)J e -^Mt)J" e -J T AHt)j 

x g s (tf(t)J + J r h(t), J*h(t) + h\t)J) (7.2) 
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for the same correlated initial state as before, p in (4-1), where T denotes the 
transpose matrix, and 



j ] A(t)T = {j:4 



( 



( j* 



: > h\t)j = (K(t) h\{t) 



A aa (t) A ab (t) 
\A ba (t) A bb {t) 

The details of these functions are given in Appendix B. 



' J a 
\ Jb J 

(7.3) 



This characteristic functional contains different types of terms from those in 
the previous example (5.11): the counter-rotating interaction provokes "squeez- 
ing." In van Hove's limit, however, these contributions disappear. Indeed, the 
van Hove limits of the constituent functions (B.4)-(B.6) of the characteristic 
functional (7.2) are 



hm Aa(r/A 2 ) = [1 + N(u s )](l - e~ r ^ T ), ( 7 .4a) 

A — 

oo 

hm 4A bb (r/X 2 )J b = [du J 6 yi + N(u)]J btU! , (7.4b) 

A — >\J J 



lim/^T/AV^^ 2 = e -r(^)r/2 e -i[A( WS )-A( W s)]r > (7 ^ 

while all other limits vanish (see Appendix B.3), and one ends up with the 
same dynamics as the previous one (6.3) except for the frequency shift; A(a>s) 
must be substituted with A(u;s) — A(u;s), where A (a;) is defined in (B.10). The 
present example again supports the validity of the theorem proved in Article 
I [5]: (a) the initial correlation disappears, (b) the state of the total system is 
factorized at all times, and (c) the reservoir remains in the mixing state, in 
van Hove's limit. The effect of the counter-rotating interaction manifests itself 
only in the frequency shift; no other differences in the resultant dynamics from 
the previous example with the rotating- wave interaction [9]. 



Furthermore, the timescales of the factorization and of the mixing are gov- 
erned by the functions K Jb ^^{t) and $j b j b (t,t), respectively (see Appendix 
B.2); they are the same as those in the previous example [Eqs. (6.11) and 
(6.12)]. This also supports the general conclusion that the free evolution of 
the reservoir plays an essential role for the factorization and the mixing, but 
the interaction does not. The counter-rotating interaction gives rise to no sig- 
nificant effect on the factorization or the mixing. 
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8 Concluding Remarks 



We have investigated two solvable models in the light of the general theorem 
proved in Article I [5]. In both cases, we confirmed that when the initial state of 
the quantum system and the reservoir is not factorized, a correct application 
of Nakajima-Zwanzig's projection method requires that the reference state 
of the latter be mixing. In addition, close scrutiny of the solvable models 
enabled us to focus on the relevant timescales. It turns out that an effective 
factorization of the state of the total system depends on the free dynamics of 
the reservoir (responsible for mixing) as well as on the interaction. Indeed, the 
free dynamics itself is sufficient to drive a complete factorization. Moreover, 
the timescales for mixing (that in turn govern the very applicability of the 
projection method in terms of the "reference" state of the reservoir) depend 
on the "size" of local observables of the reservoir: clearly, if one has access 
to information that is distributed over larger portion of the reservoir, one 
can in general detect finer deviations from mixing. The timescales at which 
Markovianity can be considered a good approximation depend on the structure 
of the local observables that one can measure, that is on the dimension of the 
(sub)system whose evolution one wants to describe. This conclusion, physically 
sound, is in some sense a strict consequence of the philosophy at the basis of the 
C*-algebraic approach to infinite systems (in the case at hand, the reservoir, 
whose observables one can measure). 

There are other very interesting problems that we have not analyzed and that 
are related to the general features of the evolutions when it is not permissible 
to consider a factorized initial state [12,14,16,17,18,19,20]. Among others, the 
problems related to the (complete) positivity of the evolution requires addi- 
tional investigations [21,22,23,24,25,26,27,28]. Another interesting issue would 
be to discuss the applicability of this method to more articulated (and in- 
triguing) thermodynamical situations, such as those of nonequilibrium steady 
states [29], shortly discussed in Article I [5] (see Fig. 1 in Article I). It is indeed 
possible to apply the method we propose to discuss the relaxation of a system 
driven by a reservoir at a nonequilibrium steady state and this aspect will be 
discussed elsewhere [30]. 
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A Prototypes of the van Hove Limits 

The characteristic functional Q[J a , J£, Jb, j\\ t] in (5.11) is expressed in terms 
of the functions given in (5.13). The van Hove limits of these functions fall into 
the following types: by taking the weak-coupling limit A — > keeping r = \H 
finite, one obtains 

(i) G(t)e^ s * — e -iVsW2 e -iA(a, s )T } ( A _ la ) 

oo 

(ii) K fg (t)^0, * fg (t,t)^f<kjft[l + N(u>)]g u , (A.lb) 

o 

(iii) (K fg *G*K gt ,)(t)e iu * t 

-> K fg (-iu s + + )K gf ,(-iu s + o+) e -r(<-s)r/2 e -iA( ws )r > (A lc) 

t t 

(iv) A 2 J dt'J dt"(K fg *G)(t-t')* gg (t\t")(G**K} !g )(t-t") 



- K fg (-iu s + + )[K fg (-iu s + 0+)]*[l + N(u s )](l - e- T ^), 

(A.ld) 

t 

(v) Jdt'^j bg (t,t')(G**K} g )(t-t') 



o 

oo 



where T(u) and A (a;) are defined in (6.4), and Kf g (s) is the Laplace transform 
of Kf g (t) in (5.14). Let us prove these results. 

(i) The van Hove limit of G(t), which is defined in (3.7), is the ordinary one 
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[9]: 



G(t)e^=r^- ^ ^fdl (A2) 

7 2m 5 + fr(A 2 § - iu s ) J 2m§ + K(-iu Q + 0+) 

which results in (A. la), by noting the formula for K(s) in (3.7), 

k(-iu s + + ) = ^T(lu s ) + tA(u s ), (A.3) 

with T(u) and A (a;) defined in (6.4). 

(ii) The van Hove limits in (A. lb) are just the long-time limits and are due to 
Riemann-Lebesgue's lemma. The timescales of the decays are determined by 
the band widths of their Fourier transforms. See Eqs. (5.14) and (6.12). 

(iii) In terms of the inverse Laplace transform, the convolution in (A.lc) is 
written as 

(K f9 * G * K gf ,)(t)e^ = I — K ^ S)K ^'[ S) e ("**)«, (A.4) 



C B 



whose van Hove limit proceeds like in (A. 2). 

(iv) Notice first that the contribution of N U u' to Eq. (A. Id) through the func- 
tion &f g (t,t'), which represents the effect of the local perturbation Lb for 

in (4.8), decays out in van Hove's limit, since the van Hove limit of this 
contribution is a generalization of (iii) but with a vanishing prefactor A 2 in 
(A. Id). Therefore, the main contribution comes from the mixing state through 

(t - t') defined in (6.13): by noting that 



oo 

= -j ^/^[l + JVHle-^H) 
o 

oo 

= / du> faa+^gx^+u* [1 + N(X 2 Co + uj s )] e -^-r') 



oo 



^ / duf^l + NMje-^-^ 

— oo 

= 2nf: s g LUs [l + N(uJs)]5(r-r > ), (A.5) 
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Eq. (A. Id) is deduced via 



T T 

J dr' J rfr' / (^ /9 *G)((r-rO/A 2 )e^( T - T ')/ A2 $W((r , -r ,/ )/A 2 )e i ^( T '' T '')/ A2 
o o 

x (G* * K}, g )((r - T")/\ 2 )e'^ T - T "^ x2 

T 

^ K fg (-iu s + + )[K f/g (-iu s + + )]*rVs)[l + N(u s )} J dr'e- r ^'. 

(A.6) 



(v) While the contribution of A/La;' decays out in van Hove's limit, which is 
shown by generalizing (ii) and (iii), that of A/^°l, 



oo 

J dcoJ* b Jl + N(u;)}g u 



ds Kfg(s) e (s+iu)T/A2 V 



OO / 



2ni s + iu s + X 2 K(s) s + iu 
j 'duJ* b Jl + N(u)] gtJ ( J 

(A.7) 



ds Kf g (\ 2 s — iu) e ST 

2ri A 2 s - i{u - u s ) + X 2 K(X 2 S - iu) T 



yields (A.le). 

The prototypes (i)-(v) lead to the van Hove limits of the components (6.2). 



B Solution to the Model with the Counter-Rotating Interaction 



We summarize the exact solution to the model with the counter-rotating in- 
teraction (7.2). 
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B.l Heisenberg Operators 



The exact solution to the Heisenberg equations of motion for a(t) = e lHt ae lHt 
and b u (t) = e iH %e~ iHt reads 



i(t) = [F(t) + X 2 F(t)]a + X 2 F(t)a^ + A j dt' F(t - t')[B(t') - B*(t% 

o 

(B.la) 

t 

b u (t) = - X J dt'e-^-^g^ait') + a f (0]> ( B - lb ) 



o 



where B(t) is defined in (3.5) and 

ds s — ius 



F(t) = [ S '^ S s e«, (B.2a) 

Cb 



with 



= - / ±L lL{S) , e* (B.2b) 



2w 



L(s) = - doj\ 9uj \ 2 ^—^. (B.3) 



s 2 + w 2 
u 

Note that F(0+) = 1, F(0+) = -iuj s and F(0+) = 0, F(0+) = 0. 
5.2 Characteristic Functional 



The characteristic functional of the state of the total system, p(t), is given by 
(7.2), which is composed of the functions 



A aa (t)^[l-\F(t) + X 2 F(t)\ 2 -X 4 \F(t)\ 2 } 
t t 

+ X 2 J dt' J dt" F{t - t') Re^ gg (t',t")F*(t - t"), (B.4a) 
o o 
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4A bb (t)J b = $ JbJb (t,t) - -\ 2 [\(K Jbg *F)(t)\ 2 + \(K* Jbg *F)(t)\ 2 ] 

t 

+ 2A 2 Im J dt' (K Jbg *ImF)(t-t')^ gJb (t',t) 
o 

t t 

+ 4A 4 J dt' J dt" {K Jbg * Im F) (t - t') 
o o 

xRe^ gg (t',t")(lmF*K* Jb9 )(t-t"), 

(B.4b) 



■AA ba (t) = \KK. hg * F)(t)F*(t) + i\\K Jb9 * lmF)(t)F*(t) 



+ \\jdt'^ Jbg {t,t')F\t-t') 



o 

t t 



+ 2A 3 y dt' J dt" (K Jb9 *lmF)(t - t')lm$P gg (t' ,t")F*(t - t"), 
o o 

(B.4c) 



AaW = -A 2 [F(t) + A 2 F(t)]F*(t) 

t t 

- ^A 2 J dt' J dt"F*(t - t') Re^ gg (t',t")F*(t - t"), (B.5a) 
o o 

JlA bb (t)J* b = -\\ 2 {K Jb9 * F)(t)(K Jb9 * F*)(t) 
t 

+ a 2 | dt' tfjjt, t') (Im F * K Jb9 ) (t - t') 
o 

t t 

+ 2A 4 J dt' j dt" (K Jbg * Im F) (t - t') 
o o 

x Re^Jt',t")(lmF*K Jbg )(t-t"), (B.5b) 
22 



1 1 r 

JlA b a(t) = --\ 3 (K Jbg *F)(t)F*(t) - -\J dt'^ Jbg (t,t')F(t-t') 



t t 



+ t\ 3 J df J dt"(K Jbg *lmF)(t-t')^ gg (t',t")F(t-t"), 
o o 

(B.5c) 



= lx[F(t) + \ 2 F(t)](K Jb9 * F*)(t) 
t t 

+ t\ 3 J dt> J df Fit -t')^ gg (t',t")(K Jbg *lmF)(t -t"), 

(B.5d) 



o o 



and 



h a (t) = F(t) + X(F * K g{m )(t) + \ 2 F(t), (B.6a) 
4h(t) =K Jbm) (t)-\(Kj bg *F)(t)-2t\ 2 (Kj bg *lmF*K gm) )(t), (B.6b) 



h a (t) = —X(F * K* g{m )(t) - X 2 F*(t), (B.6c) 
4h(t) = -\(K Jbg * F*)(t) + 2i\ 2 {K Jb9 * ImF * K* g{m )(t) : (B.6d) 

where 



oo oo 



tf g (t, t') = JdooJ dJ /:e-^(l + 2M) uuf e^g u ,. (B.7) 







5.5 Van Hove's Limit 



In addition to the prototypes (A.l), the following limits are necessary for 
the van Hove limit of the characteristic functional (7.2): by taking the weak- 
coupling limit A — > keeping r = A 2 t finite, we have 



F(t)e iu > st = [ — 
w J 2-Ki 



ds \ 2 s — 2icu s ^ T 



2m A 2 s 2 - 2iu; s s + 2u s L(\ 2 s - icu s ) 
27ri 5 + iL(-iuj s + 0+) 

(B.8a) 
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F(t)e- tuJst = I — e ST ^ 0, (B.8b) 

' 2ni + 2ilu s § + 2lu s L(\ 2 § + iu s ) 



' 2m A 2 s 2 =F 2iu s s + 2u s L(\ 2 s =F iou s ) 



where 



:b.8c) 



±iL{Tico s + + ) = ^T{u s ) ± z[A(cu s ) - A{u s )} (B.9) 



with r(a;) and A(a;) in (6.4), and 



' 2tt u + V ; 



o 



Then, the van Hove limits of the components (B.4)-(B.6) of the characteristic 
functional (7.2) yield (7.4). 



References 

[1] H. Spohn, Rev. Mod. Phys. 52 (1980) 569. 

[2] R. Kubo, M. Toda, N. Hashitsume, Statistical Physics II: Nonequilibrium 
Statistical Mechanics, 2nd ed., Springer, Berlin, 1995. 

[3] U. Weiss, Quantum Dissipative Systems, World Scientific, Singapore, 1993. 

[4] CW. Gardiner, P. Zoller, Quantum Noise, 2nd ed., Springer, Berlin, 2000. 

[5] S. Tasaki, K. Yuasa, P. Facchi, G. Kimura, H. Nakazato, I. Ohba, S. Pascazio, 



quant-ph/0602184 (2006). 

[6] S. Nakajima, Prog. Theor. Phys. 20 (1958) 948; 
R. Zwanzig, J. Chem. Phys. 33 (1960) 1338. 

[7] F. Haake, in: G. Holder (Ed.), Quantum Statistics in Optics and Solid-State 
Physics, in: Springer Tracts in Modern Physics, vol. 66, Springer, Berlin, 1973, 
pp. 98-168. 

[8] L. van Hove, Physica 21 (1955) 517; 
L. van Hove, Physica 23 (1957) 441; 
E.B. Davies, Commun. Math. Phys. 39 (1974) 91; 

E.B. Davies, Quantum Theory of Open Systems, Academic Press, London, 1976; 



24 



P.F. Palmer, J. Math. Phys. 18 (1977) 527; 

H. Spohn, J.L. Lebowitz, Adv. Chem. Phys. 38 (1979) 109; 

I. Ojima, J. Stat. Phys. 56 (1989) 203; 

L. Accardi, Y.G. Lu, I. Volovich, Quantum Theory and Its Stochastic Limit, 
Springer, Berlin, 2002. 

[9] P. Facchi, S. Pascazio, Physica A 271 (1999) 133. 

[10] G.W. Ford, J.T. Lewis, R.F. O'Connell, Phys. Rev. A 37 (1988) 4419. 

[11] G.W. Ford, R.F. O'Connell, Phys. Rev. A 61 (2000) 022110. 

[12] M. Rosenau da Costa, A.O. Caldeira, S.M. Dutra, H. Westfahl, Jr., Phys. Rev. 
A 61 (2000) 022107. 

[13] R. Haag, Local Quantum Physics: Fields, Particles, Algebras, 2nd revised and 
enlarged ed., Springer, Berlin, 1996; 

O. Bratteli, D.W. Robinson, Operator Algebras and Quantum Statistical 
Mechanics 1-2, 2nd ed., Springer, Berlin, 2002. 

[14] H. Grabert, P. Schramm, G.L. Ingold, Phys. Rep. 168 (1988) 115, and references 
therein. 

[15] G.W. Ford, M. Kac, P. Mazur, J. Math. Phys. 6 (1965) 504; 
P. Ullersma, Physica A 32 (1966) 27, 56, 74, 90; 
B.L. Hu, J.P Paz, Y. Zhang, Phys. Rev. D 45 (1992) 2843; 
B.L. Hu, A. Matacz, Phys. Rev. D 49 (1994) 6612. 

[16] G.W. Ford, J.T. Lewis, R.F. O'Connell, Phys. Rev. A 64 (2001) 032101; 
G.W. Ford, R.F. O'Connell, Phys. Rev. D 64 (2001) 105020; 
G.W. Ford, R.F. O'Connell, Ann. Phys. (N.Y.) 319 (2005) 348. 

[17] F. Haake, R. Reibold, Phys. Rev. A 32 (1985) 2462; 

A. Suarez, R. Silbey, I. Oppenheim, J. Chem. Phys. 97 (1992) 5101. 

See also, F. Haake, M. Lewenstein, Phys. Rev. A 28 (1983) 3606; 

U. Geigenmiiller, U.M. Titulaer, B.U. Felderhof, Physica A 119 (1983) 41. 

[18] V. Romero-Rochin, I. Oppenheim, Physica A 155 (1989) 52; 

V. Romero-Rochin, A. Orsky, I. Oppenheim, Physica A 156 (1989) 244. 

[19] V. Gorini, M. Verri, A. Frigerio, Physica A 161 (1989) 357. 

[20] H. Zoubi, M. Orenstien, A. Ron, Ann. Phys. (N.Y.) 313 (2004) 72. 

[21] P. Pechukas, Phys. Rev. Lett. 73 (1994) 1060; 
R. Alicki, Phys. Rev. Lett. 75 (1995) 3020; 
P. Pechukas, Phys. Rev. Lett. 75 (1995) 3021. 



25 



[22] G. Lindblad, J. Phys. A 29 (1996) 4197; 
G. Lindblad, J. Math. Phys. 39 (1998) 2763; 
A.J. van Wonderen, K. Lendi, J. Phys. A 33 (2000) 5757. 

[23] A. Royer, Phys. Rev. Lett. 77 (1996) 3272; 
A. Royer, Phys. Lett. A 315 (2003) 335. 

[24] P. Stelmachovic, V. Buzek, Phys. Rev. A 64 (2001) 062106; 
D. Salgado, J.L. Sanchez-Gomez, |quant-ph/0211164| (2002); 
P. Stelmachovic, V. Buzek, Phys. Rev. A 67 (2003) 029902(E). 

[25] H. Hayashi, G. Kimura, Y. Ota, Phys. Rev. A 67 (2003) 062109; 

D.M. Tong, J.L. Chen, L.C. Kwek, C.H. Oh, |quant-ph/0311091| (2003); 

D.M. Tong, L.C. Kwek, C.H. Oh, J.L. Chen, L. Ma, Phys. Rev. A 69 (2004) 
054102; 

D. Salgado, J.L. Sanchez-Gomez, M. Ferrero, Phys. Rev. A 70 (2004) 054102. 

[26] F. Benatti, R. Floreanini, R. Romano, J. Phys. A 35 (2002) 4955; 
F. Benatti, R. Floreanini, R. Romano, J. Phys. A 35 (2002) L551. 

[27] KM. Fonseca Romero, P. Talkner, P. Hanggi, Phys. Rev. A 69 (2004) 052109. 

[28] T.F. Jordan, A. Shaji, E.C.G. Sudarshan, Phys. Rev. A 70 (2004) 052110. 

[29] D. Ruelle, J. Stat. Phys. 98 (2000) 57; 

W. Aschbacher, V. Jaksic, Y. Pautrat, C.-A. Pillet, mp_arc 05-207 (2005); 



S. Tasaki, J. Takahashi, [cond-ma t/0606259 (2006). 
[30] S. Tasaki et al, in preparation. 



26 



